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ABSTRACT 

Magnetic field line motion and its relationship with particle 
E x B drifts are considered. It is shown that a unique and significant 
field line identification can be made whenever all field lines in a 
region pierce a conducting, but arbitrarily shaped surface. If the 
region in question is free space, particle electric drifts, in general, 
bear no relation with field line motion identified in this fashion. 
However, with a conducting plasma filling the region and allowing no 
parallel component of electric field ( E • B = 0), particles drift so 
as to always remain on the same magnetic line of force. These points 
are specifically illustrated for a time -dependent model of the earth 1 s 


magnetosphere . 



INTRODUCTION 


- 2 


A widely used concept in the electrodynamics of a guiding center plasma 
is that of frozen-in magnetic field line motion. This concept asserts 
( Newcomb, 1958) that the velocity field constructed instantaneously of 
particle E x B drift velocities is, in the appropriate circumstances, a 
legitimate (flux preserving) field line motion. 

Theoretical analyses of particle diffusion in the earth 1 s magnetosphere 
have used this concept ( Kellogg, 1959 ; Parker , i960; Davis and Chang , 1962 ; 
Nakada and Mead , 1965; Falthammar , 1965; Conrath, 1967)- During model, 
time -dependent, magnetospheric disturbances, a magnetic field line is 
identified (in these studies) by the time invariant longitude and latitude 
with which it enters the earth's dipole. Each field line is consequently 
traced through space (using the field line equations) as a function of 
time. The frozen -in field line concept is then applied in the identification 
of particle E x B drifts with the motion of field lines as so defined. 

Oar interest in this problem was stimulated by Kavanagh. 1 s (1967) recent 
criticism of the work of Nakada and Mead. Kavanagh 1 s criticism is equally 
applicable to all of the previously cited work and is based upon a direct 
calculation of particle electric drifts, using a vacuum electric field 
characteristic of the time dependent magnetospheric models. Kavanagh 
finds that the particle drifts which he calculates differ grossly from 
the field line motions of the earlier work. Seemingly, either Kavanagh 
or the earlier authors are in error. 

In this paper we show how a magnetic field line velocity of significance 
can be uniquely defined whenever all magnetic lines cut an arbitrarily shaped 
but perfectly conducting surface. It Is found that particle E x B drifts 
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coincide with the field line motion whenever no parallel component of 
electric field exists ( E • B = 0). Our results indicate that the field 
line identification and the frozen-in field assumption previously used 
are indeed valid, but in a model which differs from that assumed by 
Kavanagh by the presence of a conducting plasma. We further show 
explicitly how the electric field found by Kavanagh is modified by the 
presence of this conducting plasma, the ensuing E x B motions being then 
identical with those derived by the earlier workers. 

Field Line Identification In the Presence of a Conducting Surface 


Consider a volume of space V in which there exists a magnetic field 
B(r, t) . Let the field be sufficiently regular that it may be described 
by Euler potentials a(r, t) and g(r, t), where B _ x Vg. Define the 

rv/ 

magnetic vector potential as A _ aVg. 


The choice of a and. g is not unique, for any transformation 
a = a (a, g, t), g' - g' (a, g, t) leads to equally acceptable Euler 
potentials provided that the Jacobian, J ' " e ^~ the transformation 
is unity. Such transformations of Euler potentials are completely equiv- 
alent to gauge transformation of the vector potential A -» A + VJl • 


Now consider a stationary surface S sufficiently large that every 
line of force cuts S at least once. A choice of a and g on S (subject, 
of course, to B _ fa x Vg) fixes the gauge throughout V. 

The electric field E(r, t) is given by 
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¥e have introduced into the last form of Eq. (l) the velocity 
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It may be verified by direct substitution that v satisfies the 
equation 


v x (£ + 

and is therefore a perfectly valid flux preserving field line velocity. 
Further, with this field line identification, lines move in such a way 
as to preserve their a - p labels. This follows from the easily proven 


w x B | _ 
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relations 
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Now let us demand that the surface S be a perfect conductor so 
that we have the requirement that the normal component of B at the 
surface be constant in time. 
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Here n is the unit normal to the surface, X^ and X 2 are general coordinates 
defined on the surface, f(Xj, X 2 ) is a function which depends on the 
curvilinear nature of the X x , X 2 coordinates, and J is the Jacobian. 
Equation 6 implies J = J and hence J ( - 1 ^ = 


^ X^HO / J C = 1 ' 
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Any changes with time of a anc *- 3 on S thus correspond to relabeling 

transformations, which we eliminate by choosing ^7* _ = 0 on the 

ot ot 

surface S. 

A unique identification of magnetic field lines is thus established: 

a field line is permanently labeled by the time independent values of 

a an d g at the fixed point (or points) where it crosses S. This is, in 

fact, the identification of field lines made by the earlier authors, the 

role of the conducting surface being here assumed by the requirement that 
the same field line always enter the earth origin at the same longitude 

and latitude. 

The question now arises as to what relationship this motion has to 
the particles which might be trapped in the magnetic field. The drift 
of low energy trapped particles perpendicular to B is predominantly the 

A-' 

E x B drift velocity 
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(Magnetic gradient and curvature drifts are proportional to the energy 
of a pdrticle . ) This drift velocity is itself an acceptable field line 
velocity only if 



V y 


re b b ' 

- 


-O , 


( 8 ) 


- 7 - 


Our (now uniquely specified) field line velocity w differs from this 
drift velocity by 


v-w- 4 1 * . 


( 9 ) 


Since the tangential electric field vanishes at a conducting surface 

and ^ = 0 on S by construction, we see from Eq. (l) that (cp + \|r) is 

constant everywhere on S. Without further restriction on the electric 

B x V(cp + t|, ) 

field in the volume V , the term c — — does not, in general, 

B 2 

vanish; E is yet undetermined; and v and v differ. It is just this 
difference between v and w which is the object of Kavanagh 1 s criticism. 

Impose now the condition that E • B _ 0 by filling the volume V 
with a conducting plasma. We find from Eq. (l) that (cp + \(r) is constant 
along each field line. Furthermore, since all field lines intersect S 
and (cp + i|r) is constant on S, we conclude that when E • B = 0, (cp + if) is 
constant throughout space. From Eq. ( 9 ) it is then immediately evident 
that the electric drift velocity v and this uniquely defined field line 
velocity w are identical. The electric field is directly given as 


E 


_ i f M _ M 

- c L at at 


This assumption that E * B = 0 distinguishes 


the models developed in previous work from Kavanagh' s vacuum configuration. 

-> 

Application to a Model Magnetosphere 

In illustration of the field line identification proposed in the 
previous section, we here consider a model for time dependent magnetospheric 
distortions. The model is illustrated in Figure 1 and consists of a 
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spherical perfectly conducting earth of radius a and a dipole of moment 

M. offset at the distance x(t). Internal to the sphere it is assumed 

2ll 

that there exists the static, uniform magnetic field B._ e . The 

3 Z 

a 

conducting sphere is the surface at which field line identification is 
made . 

The magnetic field external to the sphere has two sources : the 

moving dipole, instantaneously at Z, and surface currents, both permanent 
and induced, on the sphere. At points r « j£, this magnetic field has 
the form 



Note that the undistorted magnetic field ( l co) is also dipole in nature. 
Longitudinally asymmetric compressions of the magnetosphere are caused by 
dipole motion toward the sphere ^ < cT), while expansions are the results 
of the reverse dipole motion > cT). 

It can be verified by direct differentiation (B = % x Vp) that this 
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magnetic field is represented to the indicated asymptotic order by 
the Euler potentials 



(11a) 


(lib) 


Note that at the surface of the sphere, a and p reduce to the time invariant 
values a - - — sin 2 ijj , p - 0. The unique identification of field lines 
introduced in the previous section becomes for this example a question of 
mapping in the region a < r « i, the curves <y = const., p - const, as a 
function of time (with, of course, a and p defined by Eq. 11 ) . 

Associated with the time varying magnetic field Eq. 10, is an electric 
field. Assuming the region surrounding the sphere to be a vacuum, we obtain 
an analytic representation for this field by solving the equation 
1 dB 

V x E = - — — , subject to the boundary conditions (l) that the tangential 

'w* - c at 

electric field be zero on the surface of the sphere and (2) that in the limit 



10 


of vanishing sphere radius (a - 0) the field be simply represented as 

E _ — — e x B. Such an electric field is 

- c dt x — 



( 12 ) 


Note from Eqs. (10) and (12) that at all longitudes 0, low latitude 
mirroring particles excute E x B motions away from the sun during magneto- 
spheric compression — - < 0) and toward the sun during magnetospheric 

expansion. One can directly verify from Eqs. (ll), however, that all 
magnetic field lines at low latitudes move toward the earth during the 
compressional phase. Herein is exemplified the difference between w and v 
discussed in the previous section. 
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That the electric drift is not a legitimate field line velocity 
in this case may be verified directly: 



The right hand side of Eq. (13) is, in general, non-vanishing in the region 
a < r « 

Witness now the result of adding a plasma to the region exter nal to 
the sphere. Charge separation in the plasma is assumed the source of an 
electrostatic field which cancels the component of E (as given by Eq. 12) 

(ji 0 

parallel to B. So long as — « c, currents induced in a low p plasma are 
negligible, and the magnetic field and Euler potentials as given by Eqs. 
(10) and (ll) are virtually unaltered. 

Thus, in the model as modified by this presence of plasma, the Euler 
potentials may be used to directly calculate the electric field 



- 12 - 



Note that to leading order in (~) an azimuthally symmetric lougitudinal-i./ 
directed electric field from Eq. (14) combines with the undistoried dipole 


magnetic field component of Eq. (10) to drive particles - together with 
field lines - toward the earth during magnetospheric compression. 

The electric field as given by Eq. (l J 0 is, of coarse, not, divergence- 
free. Poisson's Equation may be used to deduce the plasma charge density 
needed to effect the drastic change from the vacuum electric con figuration 
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For values representative of the earth' s magnetosphere (i = 10 26 gauss cm 3 , 

l ~ lCdo cm (~ 20 earth radii), ~ = -' - X - — — — I the maximum 

’ dt sec J 

charge density needed at six earth radii in the equatorial plane is 

approximately 3 x 10 20 . This corresponds to an order of magnitude 

(cm) 3 

difference in the number densities of protons and electrons of one part 
in 10 10 . 

Summary 

Our main point in this paper has been to indicate that in any region 
permeated by a magnetic field B, the presence of a conducting surface 
through which all field lines pass permits a unique, legitimate, and signi- 
ficant identification of the motion of magnetic field lines. The field 
line motion so defined is the E x B particle drift motion if E • B _ 0; 
otherwise the two are different. For the case E ■ B _ 0 low energy 
particles may be regarded as moving frozen to the field lines as identified 
by this prescription. 

In illustration of these points a time dependent model magnetosphere 
has been considered. Distinction between the situations where the magneto- 
sphere is free space and filled with plasma is made. The latter case is 
found to correspond to models adopted in diffusion analyses where the frozen 
in field assumption is - justifiably - used. 
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Figure Captions 

Figure 1. Geometry of the Compressed Magnetosphere. 
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